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Abstract 

Fukui et al. (2006) discovered two molecular loops in the Galactic center and ar- 
gued that the foot points of the molecular loops, two bright spots at both loops ends, 
represent the gas accumulated by the falling motion along the loops, subsequent to 
magnetic notation by the Parker instability. We have carried out sensitive CO obser- 
vations of the foot points toward /=356° at a few pc resolution in the six rotational 
transitions of CO; 12 CO(J=1-0, 3-2, 4-3, 7-6), 13 CO(J=1-0) and C 18 O(J=l-0). 
The high resolution image of 12 CO ( J=3-2) has revealed the detailed distribution of 
the high excitation gas including U shapes, the outer boundary of which shows sharp 
intensity jumps accompanying strong velocity gradients. An analysis of the multi-J 
CO transitions shows that the temperature is in a range from 30-100 K and density 
is around 10 3 -10 4 cm -3 , confirming that the foot points have high temperature and 
density although there is no prominent radiative heating source such as high mass 
stars in or around the loops. We argue that the high temperature is likely due to the 
shock heating under C-shock condition caused by the magnetic notation. We made 
a comparison of the gas distribution with theoretical numerical simulations and note 
that the U shape is consistent with numerical simulations. We also find that the region 
of highest temperature of ~ 100 K or higher inside the U shape corresponds to the 
spur having an upward flow, additionally heated up either by magnetic reconnection 
or bouncing in the interaction with the narrow neck at the bottom of the U shape. 
We note these new findings further reinforce the magnetic floatation interpretation. 

Key words: ISM: clouds — ISM: magnetic fields — magnetic loops — Radio lines: 
ISM 

1. Introduction 

The central molecular zone (hereafter CMZ, Morris & Serabyn 1996) is located in the 
inner 300 pc of the Galactic center and contains the Sgr A and Sgr B2 molecular clouds, two 
outstanding features in the CMZ (e.g., Scoville, Solomon & Penzias 1975; Fukui et al. 1977; 
Giisten & Henkel 1983). The molecular gas in the CMZ is characterized by high kinetic 
temperature from 30 K to 300 K (e.g., Rodriguez- Fernandez et al. 2001; Hiittemeister et al. 
1993; Martin et al. 2004; Oka et al. 2005; Nagai et al. 2007), and high density around 10 4 cm" 3 
(e.g., Tsuboi, Handa & Ukita 1999). The molecular gas in the CMZ also shows violent motions 
with an rms velocity dispersion of 15-30 km s _1 , much larger than those of the typical molecular 
clouds of few km s _1 outside the central kpc (Morris & Serabyn 1996; Giisten & Philipp 2004). 
It has been suggested that supernova explosions may be responsible for these properties, but 
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the total star formation efficiency is too small to account for the required abundance of OB 
stars (Morris & Serabyn 1996). The origin of these two peculiar properties, high temperatures 
and violent motions, has been puzzling since its discovery in the 1970's. It is important to 
understand the physical properties of the molecular gas to understand star formation in the 
Galactic center and, consequently, the evolution of the Galaxy. 

There are several molecular features with very broad velocity widths outside the CMZ 
including Clumps 1 and 2 (Bania 1977) and the I = 5.5° cloud (Bitran et al. 1997). These clouds 
are located outside the CMZ, perhaps distributed in the central 1 kpc, but they have not yet 
been given as much attention as in the CMZ. 

Recently, Fukui et al. (2006; hereafter F06) discovered loop-like molecular features, 
hereafter loops 1 and 2, toward /=355°-359° outside the CMZ (Figure 1) and showed that the 
two loops each have two foot points on the both ends; the foot points are bright molecular 
condensations and show velocity widths as large as ~ 50 km s -1 . F06 suggested that the 
two loops are formed by magnetic flotation driven by the Parker instability and showed that 
numerical magneto-hydrodynamics (MHD) calculations reproduce the loops successfully. This 
discovery is the first observational verification of the Parker's instability in the Galactic scale 
40 years after the prediction by Parker (1966). In the scenario of F06, the two foot points are 
interpreted as the accumulated gas formed by the down-falling motion by the stellar gravity. 
Details of such a signature has already been shown by numerical simulations by Matsumoto et 
al. (1988). F06 argues that the relatively strong magnetic field of 150 fiG, which is eventually 
a result of the strong gravitation in the inner 1 kpc of the Galaxy, makes it possible to create 
the loops throughout the Galactic center. 

Subsequently, Torii et al. (2009) made a detailed analysis of the NANTEN 12 CO( J=l- 
0) and 13 CO(J=1-0) datasets and revealed further details including helical distributions in 
loops 1 and 2 as well as associated Hi and dust features. Fujishita et al. (2009) presented the 
discovery of loop 3, another magnetically floated loop, in the same direction with loops 1 and 2 
in a positive velocity range. Machida et al. (2009) presented a three-dimensional global MHD 
simulations of the central 2-kpc magnetized gas disk, and Takahashi et al. (2009) detailed two- 
dimensional local MHD simulations of the magnetic loops. These follow-up works offer further 
supports for the magnetic flotation picture. 

In the magnetic flotation picture, the floated gas falls back down to the galactic plane 
along the loop driven by stellar gravity, and then collides with the nuclear gas disk. If the 
falling speed exceeds the sound speed, shocks must occur, leading to violent gas motions with 
enhanced density and kinetic temperature caused by compression and heating by the shocks. 
This is a conversion of the magnetic energy into the kinetic and thermal energy of the gas, and 
the magnetic flotation may provide an interpretation of physical conditions and kinematics of 
the central molecular gas. In order to test this scenario, it is important to reveal the distribution 
of temperature and density of the loop, in particular, in the foot points. 
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In this paper, we present and discuss the results of pc-scale observations in six rotational 
transitions of 12 CO (J=l-0, 3-2, 4-3, 7-6), 13 CO(J=1-0) and C 18 O(J=l-0) toward two of the 
foot points with the ASTE, Mopra, NANTEN2 mm/sub-mm telescopes. We adopt a distance 
of 8.5 kpc to the Galactic center. Details of the observations are given in section 2 and the 
results are given in section 3. An analysis of the line radiative transfer is given in section 4, a 
discussion is in section 5, and the conclusions are presented in section 6. 

2. Observations 

We observed 12 CO(J=1-0), 13 CO(J=1-0) and C 18 O(J=l-0) emission lines with the 
Mopra 22 m telescope and the 12 CO(J=3-2) line with the ASTE 10m telescope. In addition 
we observed the 12 CO(J=4-3, 7-6) lines with the NANTEN2 4m telescope. The 12 CO(J=3- 
2) observing areas were selected using the NANTEN 12 CO(J=1-0) dataset (Figure 1). We 
produced a complete view of the foot points of the loops in 12 CO( J=3-2) at a 40" grid with a 22" 
beam size. The Mopra observations were limited to a subset of the area obtained with ASTE, 
but covered the brightest features associated with the foot points. NANTEN2 observations 
were carried out toward four small regions centered on J=3-2 features identified in the ASTE 
results (Figure 2). The details of these observations are as follows (see also Tables 1, 2 and 3). 

2.1. Mopra observations 

Observations of 12 CO(J=1-0), 13 CO(J=1-0), and C 18 O(J=l-0) were carried out by 
using the 22m ATNF Mopra mm telescope in Australia, during September 2007 and August 
2008. On-the-fly (OTF) mode was used, with a unit field of 4' x 4' and we scanned in both lon- 
gitudinal and latitudinal directions separately to minimize scanning effects. The telescope had 
a half-power beamwidth (HPBW) of 33" at 100 GHz, which corresponds to 1.4 pc. The typical 
system noise temperature, T sys , was 500 K in the single side band (SSB). The Mopra telescope 
was also equipped with backend system "MOPS", providing 4096 channels across 137.5 MHz in 
each of the two orthogonal polarizations. The effective velocity resolution was 0.088 km s _1 and 
the velocity coverage was 360 km s _1 at 115 GHz. MOPS enabled simultaneous observations 
of 12 CO(J=1-0), 13 CO(J=1-0), and C 18 O(J=l-0). The pointing accuracy was checked every 
1 hour to keep within 7" by observations of the 86 HGz SiO maser from AH Sco [R.A. (2000) 
= 17 h ll m 17fl6, Dec. (2000) = -32°19'30"72]. The calibration from original output data onto 
T* scale and the baseline fitting was done with the livedata task in AIPS++. The spectra 
were gridded to a 15" spacing, and smoothed with a 36" HPBW Gaussian function. We also 
smooth the channels in velocity to a 0.86 km s _1 interval for 12 CO(J=1-0) and 13 CO(J=1-0) 
and a 2.0 km s" 1 for C 18 O(J=l-0). The spectra were converted into a T mb scale by dividing 
by "extended beam" efficiency 0.55 (Ladd et al. 2005). For comparison of CO(J=1-0) with 
other excitation lines, the spectra were spatially smoothed with a 38" Gaussian function beam, 
equivalent to the beam size of NANTEN2 at 460 GHz. We finally achieved the rms noise levels 
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in the 12 CO(J=1-0), "00(7=1-0), and C 18 O(J=l-0) spectra of 0.13 K, 0.06 K, and 0.03 K, 
respectively (see also Table 1). 

2.2. ASTE observations 

Observations of the 12 CO(J=3-2) line were made using the ASTE (Atacama 
Submillimeter Telescope Experiment) 10m sub-mm telescope of NAOJ at Pampa La Bola at 
an altitude of 4800 m in Chile (Kohno 2005; Ezawa et al. 2004, 2008) for 7 days in August 
and September 2006. The HPBW was 22" at 345 GHz. We used a position switching mode 
to cover the region shown in Figure 1, in 1399 points, with a 40" grid spacing. ASTE was 
equipped with 345 GHz SIS receiver SO 345 providing a typical system noise temperature of 
190-300 K at 345 GHz for an elevation angle of 30-80 degrees in the double-side band (DSB). 
For checking the system stability and the absolute intensity calibration, M17SW [R.A. (1950) 
=18 h 17 m 30?0, Dec. (1950) = — 16°13'6"0] was observed every 2-3 hours, the absolute temper- 
ature of which was assumed to be 69. 6K (Wang et al. 1994), and we adopt a beam efficiency of 
0.6 at 345 GHz. The spectrometer comprises four digital back-end systems (autocorrelators) 
with 2048 channels. The total frequency bandwidth is 512 MHz, corresponding to a velocity 
coverage of 450 km s" 1 with a velocity resolution of 0.43 km s _1 at 345GHz. We smoothed 
the data in velocity to a 0.86 km s _1 resolution to improve the noise level and smoothed in 
space to a 38" beam size for comparison with other excitation lines. The telescope pointing 
was measured to be accurate to within 2" by radio observations of Jupiter and W Aql [R.A. 
(2000) = 19 h 15 m 23?21, Dec. (2000) = -7°2'49"8]. Finally we got an RMS noise temperature 
level of ~0.33 K in T mb . 

2.3. NANTEN2 observations 

We used the NANTEN2 4m sub-mm telescope at 4800 m altitude at Pampa La Bola in 
Chile to observe the four regions in the foot points (Figure 2). These regions are chosen based 
on the results of 12 CO(J=3-2) (see Section 3). Observations of the 12 CO(J=4-3) transition at 
460 GHz and the 12 CO(J=7-6) transition at 810 GHz were made in June 2006 and December 
2007. The HPBW in 460 GHz and 810 GHz were measured to be 38" and 22", respectively. 
The telescope was equipped with a dual-channel 460/810 GHz receiver. DSB receiver noise 
temperatures were ~250 K at 460 GHz and ~750 K at 810 GHz. The spectrometer was an 
acousto-optical spectrometer (AOS) with a bandwidth of 1 GHz of 2048 channels and the 
velocity resolution was 0.37 km s _1 at 460 GHz and 0.21 km s" 1 at 806 GHz, respectively. We 
used the OTF mode with equatorial (J2000.0) coordinates and scanned in both RA and Dec 
directions for a 2' x 2' unit field. Main beam efficiencies at 460 GHz and 810 GHz were 0.5 
and 0.45, respectively. The raw data were calibrated into Tj£ scale and then converted to T m b 
scale by dividing by the main beam efficiency. The typical rms noise fluctuations were ~0.45 
K at 460 GHz and ~0.8 K at 810 GHz at a velocity resolution of 0.86 km s" 1 . For comparison 
of 12 CO(J=7-6) with other lines, the spectra were smoothed to 38". 
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3. Results 



3.1. Distribution of the molecular gas 

Figure 2 shows the integrated intensity distribution that includes the two foot points 
of loops 1 and 2 (Figure 2a) and part of the top of loop 1 (Figure 2b) in 12 CO(J=3-2) at a 
22" ~0.9 pc resolution. The integration range in velocity, (with respect to the local standard 
of rest, LSR), is from —180 km s _1 to —40 km s -1 . The foot points are generally elongated 
vertically to the plane and the most prominent feature ranges from / = 356.15° to 356.25° and 
b = 0.8° to 1.0°. The rest of the distribution is extended with a few local maxima, e.g., toward 
(Z, 6)=(356.17°, 0.75°), (356.13°, 0.78°) and (356.25°, 1.10°). We have chosen four peaks in the 
foot points, A-D, and a peak in the loop top as listed in Table 3 for a multi-J transition analysis 
in section 4. We note the emission is stronger in the eastern half of the distribution and that 
the eastern boundary shows a sharp intensity decrease perhaps barely resolved with the present 
resolution. Figures 3a and 4a show the distribution of 12 CO(J=1-0) and 13 CO(J=1-0) in the 
same velocity range at a 2 pc resolution with a spatial coverage limited to the bright part of 
12 CO(J=3-2) and these distributions are generally similar to the 12 CO(J=3-2) distribution. 
Figures 3b and 4b show the different velocity range, —70 to —10 km s -1 , which exhibit a "U 
shape" (see section 3.2). 

3.2. Velocity distribution 

Figures 5-7 show the 12 CO(J=3-2) velocity distribution of the foot points in two ways, 
i.e., position- velocity diagrams and velocity-channel distributions every 10 km s _1 , and Figures 
8 and 9 show the labels and auxiliary lines superposed on the velocity distributions in order to 
show the components discussed below. 

Figures 5a and 5b are latitude-velocity and longitude-velocity distributions of 
12 CO(J=3-2), respectively. The main part of the foot points, hereafter "main component", is 
distributed in a velocity range from v= —150 km s^ 1 to —40 km s _1 (Figure 8). The nature of 
the weaker emission for velocities greater than —40 km s -1 was not yet discussed in F06, while 
the very narrow velocity feature at ~ 10 km s -1 is likely foreground outside the Galactic center. 
Examination of line intensity ratios like that of 12 CO(J=3-2) to ( J=l-0) indicates that most 
of the emission seen in Figure 5 probably located in the Galactic center, as suggested by their 
high excitation (see section 4). We call the low velocity feature at —30 km s _1 to km s^ 1 the 
"subcomponent" hereafter. 

We note that the subcomponent is probably linked to the main component of the foot 
points, as suggested by the connecting broad emission at b = 0.8° and by a few additional broad 
emission components at b = 0.88°, 0.94° and 1.0° (Figures 5 and 8). The subcomponent shows 
a velocity gradient from b = 0.80° to 0.90° in the opposite sense to that of the main component 
and forms a "U shape" in the v-b diagram with the main component and the broad emission 
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at b = 0.80°. Figure 6 shows the six latitude- velocity diagrams of the foot points every 0.03° in 
/. The main component and the subcomponent are clearly seen and appear linked at b = 1.0° 
(a), b = 0.9° (d,e) and b = 0.8° (d,e,f). Below b = 0.8° we also find another U shape in the v-b 
diagram and name it U shape 2 which is distributed in b = 0.7°-0.8° and v = —70 — 10 km 
s _1 (Figures 5 and 8). A protrusion toward the galactic plane is seen in I = 356.15°-356.20° 
and b = 0.7°-0.8° at v = —70 — 50 km s^ 1 . In summary the distribution of the molecular gas 
consists of two U shapes and a protrusion as indicated in Figures 7 and 9. 

Figure 7 shows the velocity channel distributions of the main components in 12 CO( J=3- 
2). The main component appears in every panel from —120 to —40 km s _1 and subcomponent 
from —30 km s" 1 to km s _1 . We recognize the intense part delineates a U shape in the bottom 
of the foot points as seen in panels from —70 km s _1 to —10 km s" 1 (see also Figure 9). This 
U shape is also seen in 12 CO(J=1-0) (Figures 3b and 4b around 6=0.8°). 

Peaks A-D are seen over a broad velocity range. Peak A is seen in two panels in Figure 
7 from —110 to —90 km s _1 , peak B in four panels from —80 to —40 km s _1 , peak C in four 
panels from —100 to —60 km s _1 , peak D in four panels from —80 to —40 km s _1 . The most 
intense region in the foot point is peak C, which is connected to the bottom of the U shape, 
and peak D corresponds to the protrusion and shows a sharp inverse-triangle like shape. Peak 
D also seems to be connected to the component at I =356.06° and b =0.76° at v — —10 km s -1 . 
These configure the another U shape. Peak B corresponds to one of the broad features shown 
in Figures 5 and 8. Other two broad features can be seen at / = 356.25° and b = 0.90° to 1.10° 
at velocity from -60 to -20 km s" 1 and at I = 356.08° and b = 0.76° to 0.80° (Figure 9). These 
two features were not observed except in 12 CO(J=3-2) emission. A general trend in Figure 7 
is that the peaks A, C and D show sharp intensity decreases toward the east and south with 
weaker "tails" toward the north. Figure 10 shows the line profile of peak C where the 12 CO 
peak velocity changes by ~ 10 km s _1 every ~20 pc with a decrease of the peak intensity by a 
factor of 4 toward the plane. Such a trend is not clearly seen for peak B, embedded in spatially 
extended emission in the panels from —70 to —60 km s^ 1 . 

3.3. Spectra presentation and line intensity ratios 

We show the line spectra obtained at the four peaks A-D as well as at the peak position 
at the top of the loop in Figure 11. All the spectra are smoothed to a 38" Gaussian beam for 
comparison with the 12 CO(J=4-3) profiles and are smoothed to 0.86 km s -1 resolution. Only 
the C 18 O(J=l-0) spectra were smoothed to 2.0 km s 1 resolution to obtain better signal-noise 
ratios. The 12 CO(J=1-0, 3-2, 4-3) and 13 CO (J=l-0) transitions are clearly detected while 
the 12 CO ( J=7-6) transition is not detected, with an upper limit of 0.76 K. The line widths of 
all the spectra are broad with velocity extents of more than 20-40 km s _1 . 

We investigate the intensity ratios, Rj/i, in order to reveal the excitation condition of 
the gas. Here Rju stands for an intensity ratio of emissions from lien % to line j. 12 CO(J=l- 
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0), 12 CO(J=3-2), 12 CO(J=4-3), 12 CO(J=7-6), 13 CO(J=1-0) and C 18 O(J=l-0) for i and j 
are represented as 1-0, 3-2, 4-3, 7-6, 13 and 18, respectively. We show three histograms of 
line intensity ratios in the four velocity ranges of —140 — 40 km s _1 (white), —40 — 20 km s~ x 
(red), —20-0 km s _1 (orange) and 0-20 km s~ x (green) and the combinations of these velocity 
ranges (black) in Figure 12. Figure 12a shows -R3-2/1-0; this indicates that the loop emission 
is characterized by high ratios peaked at 0.6 with a range from 0.3 to 0.95 at a 90 % peak 
level, somewhat smaller than i?3_2/i_o of 0.9 in CMZ estimated by Oka et al. (2007). The 
local emission has a narrow linewidths mostly in a range from to 20 km s _1 is peaked at 0.2 
with a 90 % range of 0.1 to 0.6. The other two show peaks near 0.6 while the —20-0 km s _1 
component shows another peak at 0.3 similar to the local clouds. The other two histograms, 
-Ri-o/i3) (Figure 12b) and _R 3 _ 2 /i3 (Figure 12c), also indicate that the loop component and the 
components in a range from —40 to km s -1 show ratios distinct from the local component 
from to 20 km s -1 . Figure 13 shows a histogram of the 12 CO optical depth estimated by 
taking ratio i?i_o/i3 with assumption of the abundance ratio [ 12 C]/[ 13 C]. Riquelme et al. (2010) 
estimated [ 12 C]/[ 13 C] in the loops as ~50-70. This is a typical one in the inner Galaxy and is 
higher than 24 estimated in the CMZ (i.e. Langer & Penzias 1990). Here we shall assume the 
abundance ratio [ 12 C]/[ 13 C] ~ 53 estimated by Wilson & Rood (1994). Fi gure 13 indicates 
that the gas in the loops tends to show smaller optical depths peaked at around 3-4, while the 
local component of —20 km s _1 is peaked at around 6-10. This estimate is similar to the optical 
depth derived in the CMZ by Oka et al. (1998). We summarize that the loop component is best 
characterized by high excitation conditions as indicated by the high ratio of the 12 CO(J=3-2) 
to 12 CO(J=1-0) line intensities and the smaller 12 CO optical depth. 

Figure 14 show the distributions of the ratio i? 3 „ 2 /i-o i n velocity channel distributions. 
Generally speaking, the ratio is correlated with the 12 CO(J=3-2) intensity and the ratio be- 
comes higher than 0.7 when the 12 CO(J=3-2) intensity is higher than ~44 K km s _1 . The 
most notable enhanced ratio, corresponding to peak B (broad emission), is found in Figure 
14 at (l,b) ~ (356.18°, 0.92°) for velocity range of —80 — 10 km s -1 , where the highest ratio 
is around 2.0. The secondary enhancement is seen at (l,b) ~ (356.18°, 0.7°-0. 8°) and velocity 
range of —80 — 60 km s _1 , showing a ratio around 1.0 and corresponding to peak D. 

Figure 15 shows v-b diagrams of the R^^/i-o an d is generally consistent with Figure 
14. The most enhanced component around peak B shown in Figure 14 is found toward the 
broad feature at b ~ 0.9°-0.95° for longitude range of 356.17°-356.20°. In addition, a high 
ratio feature with a very narrow velocity width of 1-2 km s" 1 is seen at v ~ —80 km s _1 and 
I = 356.2° with a weak velocity gradient from b = 0.85° to 1.05° and another similar feature is 
seen at v ~ —30 km s _1 and / = 356.15° from b = 0.80° to 0.95°. These narrow features are 
smeared out in Figure 14. 
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4. Data analysis 



4-1. Physical properties; Clumps 

We identified clumps in the following way in the 12 CO(J=3-2) integrated intensity 
distributions shown in Figure 2; 1) Find local peaks that have the integrated intensity stronger 
than half of that in peak C, because peak C shows the maximum integrated intensity level 
in the foot point. 2) Draw a contour at two thirds of the peak integrated intensity level and 
identify it as a clump unless it contains other local peaks, 3) If there are other peaks inside 
the contour, draw another contour with interval of 10 a and find a contour which has no other 
peaks inside, and 4) If we find multiple peaks with an isolated contour when we perform the 
operation (3), we should identify these all peaks as local peaks and repeat operations (2) and 
(3). 

In this way, we identified six clumps in the foot point. Then, we derived physical 
properties of the clumps. The results are summarized in Table 5. The radii of the clumps, 
r, are ~3-6 pc, and the intensity weighted velocity dispersions (FWHM), AV, at the peak 
positions are ~20-36 km s _1 . We estimate the clump masses in two ways; the mass estimated 
from r and AV by assuming virial equilibrium, M vir , and the mass estimated from 13 CO by 
assuming the local thermodynamic equilibrium (LTE). M v ; r is calculated as follows by assuming 
uniform density distribution; 

M ™ = 209 fe)(k^) 2 M ° (1) 

13 CO column densities are calculated by assuming the LTE condition to estimate the molecular 
mass. The optical depth of 13 CO, r( 13 CO), was calculated using following equation; 

T£( 13 CO) 



r( i3 CO) = -ln 



(2) 



5.29 x (J(T ex )- 0.164) 

where, T^( 13 CO) and T ex are the radiation temperature and the excitation temperature of 
13 CO, respectively. J(T) is defined as J(T) = 1/ [exp(5.29/T) - 1]. iV( 13 CO) was estimated 
from: 

W)= ^xlO» T lgg^ g (a.-). (3) 

In this study, we assume uniform T ex of 40 K as discussed in the next subsection. We assume 
two abundance ratios [H 2 ]/[ 13 CO] = 10 6 , the value of Sgr B2 (Lis & Goldsmith 1989), and 
5 x 10 5 , the average value of local clouds (Dickman 1978), to convert iV( 13 CO) into 7V(H 2 ). 
Mco of each clump is estimated by using 

M co = 2.8m H £[£ 2 ftiV(H 2 )] M Q . (4) 

while M v ; r is in the order of ~10 5 -10 6 M , M co is only in the order of 10 4 M Q . It means that 
the dynamical state of the molecular gas is considerably different from the virial equilibrium. 

9 



If we assume the uniform velocity dispersion of 20 km s 1 as a lower limit for the turbulent 
velocity, the kinetic energy of each clumps is estimated to ~ 10 50 erg. 

4-2. Physical properties; Temperature and Density 

We applied the large velocity gradient (LVG) analysis (Goldreich & Kwan 1974; Scoville 
& Solomon 1974) to estimate the physical parameters of the molecular gas toward the loop foot 
points and the loop top by adopting a spherically symmetric uniform model having a radial 
velocity gradient dv/dr. We calculate level populations of 12 CO, 13 CO and C 18 molecular 
rotational states and line intensities under these assumptions. The LVG model requires three 
independent parameters to calculate emission line intensities; kinetic temperature, Tk, density 
of molecular hydrogen, ra(H 2 ), and X {CO) /{dv/dr). X {CO) / {dv / dr) is the abundance ratio of 
CO to H 2 divided by the velocity gradient in the cloud. We use the abundance ratios [ 12 C]/[ 13 C] 
~ 53 and [ 16 0]/[ 18 0] ~ 327 (Wilson & Rood 1994) and the molecular abundance [ 12 CO]/[H 2 ] 
~ 10 -4 as a typical value for the inner Galaxy (i.e. Frerking, Langer & Wilson 1982; Leung, 
Herbst & Huebner 1984; Blake et al. 1987). We estimate the mean velocity gradient within the 
foot point as ~9.0 km s^ 1 pc" 1 (Table 5) and X (CO) /{dv/dr) to be 1.1 x 10~ 5 for 12 CO. 

In order to solve for the temperatures and densities which reproduce the observed line 
intensity ratio, we calculate chi-square defined as below; 



N-l N 



x 2 =E E 

i=l j=i+l . 



{Robs{i,j) -#lvgQU)} 



(5) 



where N is the number of transitions of the observed molecule, i and j refer to different 
molecular transitions, R Q bs{hj) is the observed line intensity ratio from transition i to transition 
j and Rhvc{hj) is the ratio between transitions % and j estimated from the LVG calculations. 
The standard deviation a{i,j) for R b s {hj) is estimated by considering the noise level of the 
observations and the calibration error. We assume that the error of calibration from T* to 
T m b is uniformly 10 % for all observations. Because we use the same calibration factor for the 
spectra obtained with Mopra, the relative calibration error has no effect when we take a ratio 
between intensities obtained with Mopra; i.e., 12 CO(J=1-0), 13 CO(J=1-0) and C 18 O(J=l-0). 
In order to reduce the error, we use the average line intensities over 10 km s~ x bins except 
for C 18 0, because C 18 has already been smoothed to be 2.0 km s _1 velocity resolution as 
discussed in section 3.3 and Figure 11, further smoothing did not give a better noise level. The 
degree of freedom, u, is defined as follows; 

V 
2 

The data used here are derived from the line profiles in Figure 11 and the ratios are estimated 
for peaks A-D and the loop top with all transitions. Therefore, the degree of freedom is given 
as 15 for N= 6. The intensities of all transitions used here are listed in Table 4. Because we 
could not detect any emissions of 12 CO(J=7-6), we use the la noise levels for calculations as 
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(6) 



upper limit. Since y 2 is weighted by a 2 , this assumption does not affect the results. 

Figure 16 shows the loci of constant Rj/i for the five Rj/i's as a function of density and 
temperature. We find #3-2/1-0 is sensitive to a large density range from 10 2 cm" 3 to 10 5 cm" 3 
while i?4_3/i_o is sensitive to higher density above 10 3 cm -3 . We also note that # 1 _ /i3 is 
nearly orthogonal to #3-2/1-0 f° r density lower than 10 4 cm -3 , making the combination useful 
in obtaining solutions. 

Figures 17a-17e and Table 6 show the results of fitting the data obtained with a x 2 
minimization approach to find the best solution of temperature and density. Each thick locus 
surrounding the cross indicates the x 2 °f 25.0, which corresponds to the 5% confidence level of 
X 2 distribution with 15 degree of freedom. The crosses denote the lowest point of x 2 ■ Three 
intensity ratios, #3-2/1-0, #4-3/1-0 an d #1-0/13, are also shown by thin lines. The —75 km 
s" 1 and —15 km s" 1 components at peak B do not have any solutions with 5% confidence 
level. Another presentation of the results is found in Figure 18 for the five peaks along with 
the 12 CO(J=3-2) line profiles. An error bar in the figure is defined as 5% confidence level 
of x 2 distribution. The density is in the range from 10 3 to 10 4 cm -3 and the temperature is 
from ~ 30-40 K to 100 K or higher. As expected from the #3-2/1-0 values in Figures 14 and 
15, the broad emission connecting different sides of the U shape (peak B) shows the highest 
temperature, 100 K or higher, amongst the five peaks. 

5. Discussion 

There are no candidates for protostars identified by IRAS in the observed regions. The 
10 GHz radio continuum emission observed by Handa et al. (1987) gives an upper limit of 0.1 
K. This 10 GHz flux density is equivalent to 10 43 ' 5 FUV photons s _1 , assuming 220 square arc 
minutes as the area of the foot point and 8000 K electron temperature. That corresponds to 
a single B3 star or a star of later spectral type (Kurtz, Churchwell & Wood 1994; Panagia 
1973; Mezger & Henderson 1967). Only one ultra compact Hll region (UCHll) is identified at 
(/, b) ~ (356.25°, 0.7°) at a velocity of +120 km s _1 from Ha recombination line observations 
(Caswell & Haynes 1987; Lockman 1989). So, the UCHll is not associated with the foot point. 
Here, there is no indication of significant formation of massive stars. Some other mechanism 
for heating the molecular gas is required, in particular, towards peaks B and C. 

The magnetic flotation model presented by F06 is able to offer another explanation to 
heat the gas by shocks. F06 suggested that loops 1 and 2 are created by magnetic buoyancy 
driven by the Parker instability. The floated gas falls down to the Galactic plane along the 
magnetic loop, and the accumulated gas forms a massive cloud in the foot point of the loop. 
In this case, if the fall down velocity of the gas exceeds the speed of sound, shock fronts are 
formed in the foot point (Matsumoto et al. 1988). The magnetic field is estimated to be 150 
fj,G in the loops under the assumption of energy equipartition, and the velocity of the gas is 
estimated to be similar to the Alfven speed, ~ 24 km s" 1 (F06). 
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The structure of shock waves in molecular clouds is an important subject. Draine, 
Roberge, & Dalgarno (1983) made such calculations including the effects of ion- neutral stream- 
ing driven by the magnetic field. They found that shock waves in molecular clouds are usually 
C-type shock waves, mediated by the dissipation accompanying ion-neutral streaming, and in 
which all of the hydrodynamic variables are continuous. In the foot point of the loops, C-type 
shock seems to be a viable model because they occur with a strong magnetic field and at a 
moderate shock speed. The limiting shock speed for C-type shocks not dissociating the H2 
molecules is estimated to be 45 km s~ x (Draine, Roberge, & Dalgarno 1983). The Alfven 
speed 24 km s _1 estimated by F06 and the velocity dispersion of the clumps shown in Table 5 
are both within this range and hence we adopt C-type shock as the heating mechanism for the 
loops. 

Draine, Roberge, & Dalgarno (1983) show that the neutral gas temperature is in a range 
100-1000 K for a density of 10 4 cm" 3 , B of 100 fiG and shock speed of 20-40 km s -1 . These 
parameters are consistent with density of 10 3 -10 4 cm -3 , magnetic field of 150 /xG and a falling 
speed of 24 km s _1 (F06). Such high temperature is consistent with the present estimate of 
temperature higher than 30-40 K, since the high temperature layer should be thin like 10 16 -10 17 
cm, ~10-100 times smaller than the current resolution of a few pc. Chieze, Pineau des Forets 
& Flower (1998) also presented results on shock parameters consistent with the above. 

We have compared the present results with local two-dimensional MHD numerical sim- 
ulations (Matsumoto et al. 1988; Takahashi et al. 2009). The distribution of the foot point 
is basically characterized by a U shape both in space (section 3.2 and Figures 3b, 4b and 7) 
and in velocity (section 3.2 and Figures 5 and 6). MHD simulations indeed indicate that the 
foot point between two loops becomes U shaped, consisting of the two down flows from the two 
loops that merges at the bottom (see e.g., Figures 4 in Takahashi et al. 2009). In velocity space, 
the foot point foreground and background components have a large velocity splitting of order 
the Alfven speed. It is preferable to compare these observational properties with the MHD 
calculations to confirm the scenario. There is, however, no previous study which focuses on the 
details in the foot point of the galactic loops. In the actual formation of loops, we need to take 
into account a time lag in falling motion between the two sides of a foot point as well as the 
other physical differences between them. In the MHD calculations, such conditions are simply 
assumed to be uniform. It is neither shown yet by MHD calculations how the compact broad 
features connecting the U shape are formed (Figures 8 and 9). In order to better understand 
the whole process, we must await detailed numerical simulations of the foot point formation 
with take into account the observed parameters and also time evolution. 

The present analysis indicates that the highest ratio in the i?3_2/i_o is formed not to- 
ward the bottom of the foot points but well above the foot points apparently separated from 
the strongly shocked layer in the bottom of the U shape. It is interesting to compare these 
results with the solar foot points for which a number of observations have been accumulated 



12 



on magnetic activity. 

Isobe, Tripathi & Archontis (2007) argued that anti-parallel field lines are produced in 
the foot point when multiple loops rose. This causes magnetic reconnection that can lead to a 
jet-like structure rising upward. Or, alternatively, in the U shaped region the field lines become 
compressed to reduce the falling gas flux, leading to gas bouncing at the narrow neck (Kudoh 
& Shibata 1999). We present a schematic figure of these scenarios in Figure 19. Numerical 
simulations in the solar case suggest that the initial condition may determine which process 
becomes dominant, either reconnection or bouncing. The basic physics is common to the spurs 
in Galactic loops as already discussed by Matsumoto et al. (1988) and Takahashi et al. (2009). 

We shall make a crude estimate of the energies concerned, both the magnetic and grav- 
itational energies released. Here we assume the concerned volume of the U shape as a cylinder 
with a radius of 10 pc and a height of 30 pc. Gravitational energy is released by the gas falling 
down along the loop that can lead to shock heating of the U shape. The heights of the tops 
of loops 1 and 2 and the U shape from the galactic plane are estimated as ~200 pc, ~300 
pc and ~130 pc, respectively, at a distance to the Galactic center of 8.5 kpc, and the total 
molecular mass in the U shape is estimated as ~ 5 x 10 4 M Q . Potential energies from the 
U shape to the tops of loops 1 and 2 divided by the time scale of 10 7 years (Machida et al. 
2009; Takahashi et al. 2009) at the galacto-centric radius of 670 pc are derived as 0.5 x 10 37 
erg s _1 and 1.3 x 10 37 erg s _1 , respectively, by using the modified Miyamoto-Nagai potential 
(Miyamoto & Nagai 1975; Sofue 1996). If the reconnection takes place, we need to take into 
account the additional magnetic energy release. We use the magnetic field strength and Alfven 
speed of 150 /xG and 24 km s _1 , respectively, estimated by F06. According to the scenario 
given above, the field lines move parallel to the galactic plane and inflow into the U shape. If 
we roughly assume the inflow speed to be 10% of the Alfven speed, the total magnetic energy 
accumulated into the cylinder is estimated to ~ 4 x 10 50 erg in ~10 7 years. This energy is 
comparable to the magnetic energy of ~ 10 51 ergs for a typical single loop given by Machida 
et al. (2009). In ~ 10 6 years after the magnetic reconnection occurs, the maximum available 
power of the reconnection is derived as ~ 1.3 x 10 37 erg s -1 . Thus, the total available energy 
injected to the U shape is 1.8 x 10 37 -2.6 x 10 37 erg s -1 by summing up both the magnetic and 
gravitational energies. The cooling power in peak B is estimated to be ~ 1.7 x 10 36 erg s _1 
for kinetic temperature of 100 K and density of 4 x 10 3 cm -3 for a uniform sphere with ~3 pc 
radius and that over the entire U shape is estimated to be ~ 4.3 x 10 36 erg s~ x for 40 K and 
4 x 10 3 cm -3 (Goldsmith & Langer 1978). These values correspond to nearly comparable to, 
or ~20% of, the maximum available power. We therefore infer that the reconnection is able to 
explain the heating of the warm gas in the U shape, while obviously we require more detailed 
model simulations and direct measurements of the magnetic field to reach a firm conclusion on 
the physical processes related to the heating mechanism. 
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6. Conclusions 



We have made detailed high-resolution observations of the foot points of the molecular 
loops 1 and 2 (Fukui et al. 2006) in six rotational transitions of the interstellar CO molecule. 
The main conclusions are summarized below; 

1) The foot points have sharp intensity gradients toward the south and east, as is con- 
sistent with the shock formation at the bottom of the foot points (Figure 2b). The foot points 
have several major peaks having 10 4 -10 5 M & gas in each. The foot points are characterized by 
two U shapes both in space and velocity and the protrusion (Figures 5-8). We suggest that the 
U shape may be formed by merging of two down flows between two loops as derived in MHD 
numerical simulations. 

2) Toward five selected peaks of 12 CO(J=3-2), peaks AD and peak top (Figures 2-4), 
we have carried out a multi-line LVG analysis of line radiation transfer and derived the molecular 
temperature and density. The four peaks in the foot point show rather high temperatures of 
~30-40 K and density of 10 3 cm -3 to 10 4 cm -3 (Figures 17 and 18). Among the four peaks, 
the peak toward the central region of the foot point, peak B, shows the highest temperature of 
~ 100 K or more. 

3) We compared the results with calculations of C-shocks by Draine, Roberge, & 
Dalgarno (1983) and find that the derived temperature and density are roughly consistent 
with theoretical estimates for magnetic field around 100 /iG, suggesting that shock heating is 
a viable explanation for the high temperature over the foot point. In addition, by comparing 
theoretical work on the solar activity, we argue that the warmest region in the central part of 
the foot point may be additionally heated up either by magnetic reconnection or by upward 
flowing gas bounced by the narrow neck in the foot point. 
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Tabic 1: Observed lines 



transition 


frequency 
(GHz) 


beamsizc 
(") 


observing mode 


noise r.m.s. 
(K) 


telescope 


12 C 16 O(J=l-0) 


115.27120 


33 


On the Fly 


0.13 


Mopra 


13 C 16 O(J=l-0) 


110.20137 


33 


On the Fly 


0.06 


Mopra 


12 C 18 O(J=l-0) 


109.78218 


33 


On the Fly 


0.03 


Mopra 


12 C 16 0(J=3-2) 


345.79599 


22 


Position Switching 


0.33 


ASTE 


12 C 16 0(J=4-3) 


461.04077 


38 


On the Fly 


0.60 


NANTEN2 


12 C 16 0(J=7-6) 


806.65181 


22 


On the Fly 


0.76 


NANTEN2 



Note. - Column (5): Typical noise r.m.s. per channel. The channel width of all the 
data is 0.86 km s — 1 . 



Table 2: Observed areas of position-switch observaions 



transition 


coord. 


observed area 


grid 
(") 


•Wpoint 


region 


12 CO(J=3-2) 


galactic 


(356.033°, 0.633°)-(356.267°, 0.811°) 
(356.267°, 0.822°)-(356.300°, 1.267°) 


40 


1399 


foot point 




galactic 


(356.462°, 1.191°)-(356.591°, 1.492) 


40 


2336 


top 



Note. - Column (2) : Coordinate system of observations. Column (4) : Grid intervals of observations. 
Column (5) : Number of total observed points. 



Table 3: Observed areas of On-thc-ny observations 









map size 
(') 


grid 






transition 


coord. 


reference position 


(") 


■*Y map 


region 


12 co, 13 co 




356.222°, 0.845° 


3 - 4 


15 


n 


foot point 


and C 18 O(J=l-0) 


galactic 


356.556°, 1.333° 


2 


15 


l 


top 






17 h 31 m 54?58, -31°08'13."76 


2 


15 


l 


top 






17 h 32 m 03?68, -31°30'55"44 


2 


15 


l 


peak A 


12 CO(J=4-3) 










12 CO(J=7-6) 


Equatorial(J2000.0) 


17 h 32 m 30?36, -31°40'04"02 


2 


15 


l 


peak B 






17 h 32 m 58?84, -31°40'59"49 


2 


10 


i 


peak C 






17 h 33 m 22?07, -31°47'51"85 


2 


15 


i 


peak D 



Note. - Column (2) : Coordinate system of observations. Column (4) : Size of the unit square of OTF scan. 
Column (5) : Grid intervals of the out put data. Column (6) : Number of maps that were observed with 
each reference position. 
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Table 4: 10 km 


s 1 average 


intensity of the peaks. 




Peak 






10 km 


1 average intensity (K) 




/l 1 \ 

(km s j 


12 CO(1-0) 


12 CO(3-2) 


12 CO(4-3) 


13 CO(1-0) 


C 18 O(l-0) 




-110 


3.31 


3.31 


2.29 


0.21 






-100 


5.26 


5.49 


4.00 


0.63 


0.05 


A 


-90 
-80 


4.27 
2.40 


3.41 
1.44 


2.24 
0.71 


0.51 
0.20 






-68 


2.06 


1.23 


0.61 


0.14 






-58 


2.43 


1.88 


0.90 


0.21 






-105 


2.39 


1.39 


0.54 


0.15 






-95 


4.75 


3.77 


1.31 


0.36 






-85 


6.56 


6.00 


3.10 


0.52 






-75 


4.52 


5.76 


3.54 


0.40 






-65 


5.24 


7.10 


4.81 


0.33 




B 


-55 


2.84 


5.48 


3.62 


0.15 






-45 


0.87 


1.94 


1.45 


0.08 






-35 


0.82 


1.51 


0.67 


0.14 






-25 


1.62 


2.50 


1.85 


0.20 






-15 


2.58 


3.42 


2.44 


0.39 






-5 


5.07 


4.30 


2.85 


0.54 






-95 


8.90 


6.21 


4.41 


0.53 


0.04 




-85 


10.33 


9.50 


7.89 


0.92 


0.08 


C 


-75 
-65 


8.22 
6.75 


7.54 
5.62 


6.05 
4.08 


0.69 
0.54 


0.06 
0.04 




-55 


5.34 


3.19 


1.87 


0.42 






-45 


4.45 


1.84 


0.94 


0.38 






-68 


6.61 


5.83 


4.02 


0.50 




D 


-58 
-48 


10.95 
5.42 


10.65 
4.75 


8.10 
3.19 


0.99 
0.47 


0.10 




-38 


2.77 


1.02 


0.77 


0.29 






-80 


3.35 


1.36 


1.00 


0.24 




TOP 


-70 
-60 


4.83 
6.23 


2.84 
4.46 


2.28 
4.06 


0.50 
0.69 






-50 


2.09 


1.66 


1.32 


0.13 





Note. " — " stands for non-detection above 3o\ Column (2) indicates Center velocity for 
calculations. 
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Table 5: Physical properties of the clumps. 





b 


AV 


r 


M vir 


M co 








(km s _1 ) 








region 


(°) 


(°) 


(PC) 


(x 10 4 M o ) 


(x 10 4 M o ) 




356.124 


0.780 


26.0 


4.6 


65.0 


1.2+ / 2.5* 




356.149 


0.913 


30.4 


2.9 


56.0 


0.7* / 1.4* 


peak B 


356.174 


0.708 


19.2 


2.8 


21.6 


0.6* / 1.2* 


peak D 


356.187 


0.988 


33.1 


2.4 


55.0 


0.4* / 0.8* 




356.220 


0.838 


36.0 


5.9 


159.8 


3.7* / 7.4* 


peak C 


356.254 


1.113 


38.2 


3.9 


118.9 


1.5* / 3.0* 


peak A 



T [H 2 ]/[ 1:i CO]= 5 X 10 5 (Dickman 1978) was used for estimation. 



* [H 2 ]/[ 13 CO]= 10 6 (Lis & Goldsmith 1989) was used for estimation. 

Note. - Column (1, 2) : Peak position of the clump. Column (3) : Intensity weighted 
velocity dispersion (FWHM). Column (4) : Radius of the clumps. Column (5) : Clump 
mass derived with assumption of Virial equilibrium. M v ; r = 209r(AV) 2 . Column (6) : 
Clump mass derived by using 13 CO(J=1-0) with LTE assumption. 
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Table 6: 


LVG results 


for X(CO)/(dv/dr) = 


1.1 x 10~ 5 . 




Peak 


v Isr 


n(H a ) (cm- 3 ) 


T k 


(TO 




(km s — M 


X 2 < 25 


min. x 2 


X 2 < 25 


min. x 2 


min. x 2 




-110 


10 3.1_ 10 4.2 


10 3.6 


38-128 


53 


2.88 




-100 


10 3.8_ 1Q 4.0 


10 3.9 


35-41 


38 


22.84 


A 


-90 


10 3.5_ 10 3.9 


1Q 3.6 


24-36 


28 


15.77 




-80 


10 2.9_ 10 3.6 


1Q 3.3 


19-65 


27 


3.24 




-68 


10 2.4_ 10 3.7 


1Q 3.2 


> 20 


32 


3.35 




-58 


10 3.0_ 10 3.8 


1Q 3.4 


22-61 


30 


4.45 




-105 


10 2.4„ 10 3.5 


10 31 


> 20 


33 


3.34 




-95 


10 3.1_ 10 3.4 


1Q 3.3 


27-40 


32 


17.16 




-85 


1Q 3.4_ 10 3.5 


1Q 3.4 


37-38 


37 


24.61 




-75 


— 


— 






30.52 




—65 


1Q 3.5_ 10 3.8 


1Q 3.7 


57-76 


65 


18.85 


B 


—55 


10 3.3_ 1Q 3.7 


1Q 3.5 


> 97 


128 


22.55 




-45 


10 3.2_ 1Q 4.3 


1Q 3.7 


> 60 


98 


19.68 




-35 


10 4.0_ 10 4.4 


io 4 - 2 


44-59 


48 


24.76 




-25 


10 3.8_ 10 4.7 


10 4.0 


41-91 


51 


20.75 




-15 


— 


— 






30.36 




-5 


10 3.5_ 1Q 3.8 


1Q 3.6 


28 - 37 


31 


19.24 




—95 


1Q3.2_ 10 3.5 


1Q 3.3 


41-56 


47 


13.13 




-85 


10 3.8_ 10 3.8 


1Q 3.8 


42-43 


42 


24.59 


C 


-75 
-65 


10 3.5_ 1Q 3.9 
10 3.4_ 10 3.8 


10 3.8 
1Q 3.6 


38-50 
34-48 


44 
40 


11.78 
11.76 




-55 


10 3.2_ 1Q 3.5 


1Q 3.3 


26-39 


31 


10.68 




-45 


10 3.1_ 1Q 3.3 


1Q 3.2 


19-28 


22 


13.85 




-68 


10 3.4_ 10 3.8 


1Q 3.5 


37-50 


42 


13.44 


D 


-58 
-48 


1Q 3.7_ 10 3.8 
10 3.4_ 10 3.9 


10 3.8 
1Q 3.6 


39-43 
32-46 


41 
37 


22.30 
13.02 




-38 


10 3.1_ 10 3.4 


1Q 3.3 


17-24 


20 


17.55 




—80 


1 2.2_ 10 3.4 


1Q 3.0 


> 18 


37 


12.88 


TOP 


-70 


10 3.2_ 10 3.8 


1Q 3.8 


20-47 


26 


10.00 




-60 


10 3.4_ 10 4.1 


1Q 3.8 


24-46 


32 


5.75 




-50 


> io 2 - 6 


1Q 3.6 


> 22 


42 


4.67 



Note. - " — " stands for no-solution above 5% confidence level of x 2 distribution. Column 
(3): Number density, n(H2) cm -3 , for that \ 2 is l ess than 25. (4) : Number density, ra(fi2) 
cm -3 , at the point that minimum \ 2 is found. (5) : Kinetic temperature, K, for that \ 2 
is less than 25. (6) : Kinetic temperature, K, at the point that minimum \ 2 is found. (7) 
: minimum \ 2 
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Fig. 1: (a-b) Integrated intensity distributions of loops 1 and 2 in 12 CO(J=1-0) obtained by 
NANTEN 4m telescope (Fukui et al. 2006). Yellow boxes in each image show the observed 
region in 12 CO(J=3-2) by ASTE. (a) Loop 1 : The integration range in velocity is from —180 
to —90 km s" 1 . Contours are illustrated from 7 K km s" 1 (white) with an interval of 50 K km 
s _1 . (b) Loop 2 : The integration range in velocity is from —90 to —40 km s _1 . Contour levels 
are the same as that in (a), (c) Schematic image of loops 1 and 2. The foot points and tops of 
the loops are depicted by circles and dashed boxes. 



22 




' ■ I ■ I ' I I > ■ ■ ■ I ■ ' ' ' I I > ■ ■ I ■ " I I 

356.3 356.2 356.1 
Galactic Longitude [degree] 



Fig. 2: Integrated intensity distributions of 12 CO(J=3-2) emissions at the top and foot point 
of the loops, (a) Top of the loopl. Integration range is from —180 to —40 km s _1 . Contours 
are plotted every 40 K km s _1 from 15 K km s _1 . (b) Foot point of the loops. Integration range 
and contour levels are the same as that in (a). The dotted box shows the observing region of 
12 CO(J=1-0) emissions and 13 CO(J=1-0) emissions, that are shown in Figures 3 and 4. The 
small boxes show the observed reigons of 12 CO(J=4-3, 7-6) emissions. 
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Fig. 3: Integrated intensity distribution of 12 CO(J=1-0) emissions at the foot point of the 
loops. Integration ranges are from —180 to —40 km s _1 (a) and from —70 to —10 km s _1 (b). 
Contours are plotted every 30 K km s _1 from 2 K km s" 1 . The labeled clumps are discussed 
in the text. 
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Fig. 4: Integrated intensity distribution of 13 CO(J=1-0) emissions at the foot point of the 
loops. Integration ranges are from —180 to —40 km s _1 (a) and from —70 to —10 km s^ 1 (b). 
Contours are plotted every 5 K km s _1 from 2 K km s -1 . 




V\sr [km s _1 ] Galactic Longitude [degree] 

Fig. 5: (a) Velocity - galactic latitude diagram of the footpoint of the loops in 12 CO(J=3-2) 
averaged from 356.03° to 356.27° in galactic longitude. Contours are plotted every 0.2 K. (b) 
Galactic longitude-velocity diagram of the footpoint of the loops in 12 CO(J=3-2) averaged 
from 0.63° to 1.27° in galactic latitude. Contours are plotted every 0.2 K. 



25 




-100 -50 -100 -50 -100 -50 

V< s , [km s- 1 ] It [km s" 1 ] Ww [km s" 1 ] 



Fig. 6: (a-f) Longitude channel maps of 12 CO(J=3-2) averaged over successive 100" intervals. 
Contours are plotted every 1 K from 0.5 K. The figure in the left side is the integrated intensity 
distributions which is the same as the image shown in figure 2. Solid lines in the figure show 
the integration ranges of figures (a)-(f). 
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Fig. 7: Velocity channel maps of 12 CO(J=3-2) integrated over successive 10 km s 1 . Contours 
are illustrated every 10 K km s _1 from 4 K km s _1 . 
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Fig. 8: Positions of main component, subcomponent and broad emission regions are superposed 
on the 12 CO(J=3-2) fi gure shown in Figure 5. 
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Fig. 9: Positions of main component, subcomponent and broad emissions are superposed on 
the velocity channel maps shown in Figure 7. Velocity range is from —70 to —10 km s" 1 . 
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Galactic Longitude [degree] V\sr [km s _1 ] 

Fig. 10: (a) Peak velocity map around peak C estimated with 12 CO(J=3-2) is shown in the 
color image. Contours show the integrated intensity levels of 12 CO(J=3-2) integrated from 
-140 to -40 km s" 1 and are plotted every 60 K km s" 1 . (b) 12 CO(J =3-2) and 13 CO(J=1-0) 
spectra of the four points on the red arrow in figure (a). The order of the spectra is the direction 
of the red arrow, from peak C to the left-bottom of the peak C. 
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Fig. 11: Observed CO spectra at the peak of each of the clumps shown in Figures 3 and 4. 
The C 18 O(J=l-0) spectra intensity scale have been multiplied by two for clearly. The vertical 
dashed lines are drawn with 20 km s _1 intervals from —120 km s" 1 . The gray scale boxes show 
the regions where the LVG calculations was carried out (Figure 18). 
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Fig. 12: Intensity weighted frequency distributions of the CO intensity ratio. The data were 
smoothed to a 60" spatial resolutions with a gaussian function and smoothed to a 2 km s" 1 
velocity resolution. The white lines show the contributions of the loops (—140 km s _1 to —30 
km s -1 ), and the red, orange and green lines show the contributions of the local components 
with a interval of 20 km s" 1 from -40 km s" 1 . (a) R 3 - 2 /i-o = 12 CO( J=3-2)/ 12 CO( J=l-0). 
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Fig. 13: Histogram of the 12 CO( J=l-0) optical depth in the foot point of the loops. The data 
are smoothed to a 5 km s _1 resolutions. The white line shows the contributions of the loops 
( — 140 km s _1 to —30 km s _1 ), and the red, orange and green lines show the contributions of 
the local components, same as in Figure 12. 
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Fig. 14: Velocity channel maps of the 12 CO(J=3-2)/ 12 CO(J=l-0) intensity ratio integrated 
over successive 10 km s _1 channels. Contours are 12 CO( J =3-2) and illustrated every 10 K km 
s _1 from 4 K km s" 1 . 
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Fig. 15: Longitude channel maps of the 12 CO(J=3-2)/ 12 CO(J=l-0) intensity ratio averaged 
for successive 45" with an interval of 120". Contours are the averaged intensity of 12 CO( J=3-2) 
and are plotted every 1 K from 0.5 K. 
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Fig. 16: Intensity ratio distributions in the density-temperature space calculated from the 
LVG model for X(CO)/(dv/dr) of 1.1 x 10~ 5 . #3-2/1-0, #4-3/1-0, #7-6/1-0, #1-0/13 and 
#i_o/i8 stand for the intensity ratio of 12 CO( J=3-2)/ 12 CO( J=l-0), 12 CO( J=4-3)/ 12 CO( J=l- 
0), 12 CO(J=7-6)/ 12 CO(J=l-0), 12 CO(J=1-0)/ 13 CO(J=1-0) and 12 CO( J=l-0)/C 18 O( J=l- 
0), respectively. 
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Fig. 17a: LVG results of peak A. Values shown at left-top of the figures show the center 
velocity of the calculated velocity range. Crosses denote the lowest point of chi-square \ 2 ■ 
Contours surrounding the cross indicate \ 2 = 25.0 which correspond to 5% confidence level of 
X 2 distribution with 15 degree of freedom. Thin lines show the typical intensity ratios; i?3_2/i_o 
(solid lines), i?4_3/i_o (dashed lines) and i?i-o/i3 (dotted lines). 
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Fig. 17b: LVG results of peak B. Details are the same in Figure 17a. 
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Fig. 17c: LVG results of peak C. Details are the same in Figure 17a. 
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Fig. 17d: LVG results of peak D. Details are the same in Figure 17a. 
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Fig. 17e: LVG results of peak top. Details are the same in Figure 17a. 
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Fig. 18: LVG results for X(CO)/(dv/dr) of 1.1 xlCT 5 for the five spectra at the peaks. The 
horizontal axis of all figures is LSR velocity. The top row shows the 12 CO(J=3-2) spectra at 
the peaks A-D, and the peak in the loop top. The second and third rows show the number 
density, n(B.2) cm -3 and the kinetic temperature, Tk K, respectively. Circles show the lowest 
point of x 2 , and the error range is defined as a 5% confidence level of x 2 distribution with 15 
degree of freedom, that corresponds to % 2 =25.0. The minimum \ 2 are shown in fourth row. 
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Fig. 19: Schematic image of reconnection at the foot point between two loops. 
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